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Abstract 

We study a system of two quantum dots connected by a hopping bridge. Both the dots and 
connecting region are assumed to be in universal crossover regimes between Gaussian Orthogonal 
and Unitary ensembles. Using a diagrammatic approach appropriate for energy separations much 
larger than the level spacing we obtain the ensemble-averaged one- and two-particle Green's func- 
tions. It turns out that the diffuson and cooperon parts of the two-particle Green's function can 
be described by separate scaling functions. We then use this information to investigate a model 
interacting system in which one dot has an attractive s-wave reduced Bardeen-Cooper-Schrieffer 
interaction, while the other is noninteracting but subject to an orbital magnetic field. We find that 
the critical temperature is nonmonotonic in the flux through the second dot in a certain regime of 
interdot coupling. Likewise, the fluctuation magnetization above the critical temperature is also 
nonmonotonic in this regime, can be either diamagnetic or paramagnetic, and can be deduced from 
the cooperon scaling function. 

PACS numbers: 73.21.La, 05.40.-a, 73.50.Jt 
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I. INTRODUCTION 



The idea of describing a physical system by a random matrix Hamiltonian to explain 
its spectral properties goes back to Wigner-i^. It was further developed by Dyson, Mehta 
and others, and became the basis for Random Matrix Theory (RMT)^. First introduced in 
nuclear physics, RMT has been used with great success in other branches of physics and 
mathematics. A notable example was a conjecture by Gorkov and Eliashberg^ that the 
single-particle spectrum of a diffusive metallic grain is controlled by RMT. This conjecture 
was proved by Altshuler and Shklovskii^ who used diagrammatic methods and by Efetov 
who used the supersymmetry method^. In 1984 Bohigas, Giannoni and Schmit^ conjectured 
that RMT could also be employed in the study of ballistic quantum systems whose dynam- 
ics is chaotic in the classical limit. Their conjecture broadened the area of applicability 
of RMT enormously and was supported by numerous ensuing experiments and numerical 
simulations^i^i^ii^. The crucial energy scale for the applicability of RMT is the Thouless en- 
ergy Et = h/Terg, where Terg is the time for a wave packet to spread over the entire system. 
For a diffusive system of size L, we have Et ^ hD/L'^, while for a ballistic/chaotic system 
we have Et — hvp/L, where vp is the Fermi velocity. 

In this paper we consider a system of two quantum dots/nanoparticles which are coupled 
by a hopping bridge. The motion of electrons inside each dot can be either ballistic or 
diffusive. In the case of ballistic dots we assume that the dots have irregular shapes leading 
to classically chaotic motion, so that RMT is applicable. 

RMT Hamiltonians fall into three main ensembles^. These are the Gaussian Orthogonal 
Ensemble (GOE), Gaussian Unitary Ensemble (GUE), and Gaussian Symplectic Ensemble 
(GSE). They are classified according to their properties time-reversal (TR). The Hamilto- 
nians invariant with respect to TR belong to the GOE. An example of GOE is a quantum 
dot which has no spin-orbit coupling and is not subject to an external magnetic field. GUE 
Hamiltonians, on the contrary, are not invariant with respect to TR and describe motion 
in an orbital magnetic field, with or without spin-orbit coupling. Hamiltonians from GSE 
group describe systems of particles with Kramers degeneracy that are TR invariant but 
have no spatial symmetries, and correspond to systems with spin-orbit coupling but with 
no orbital magnetic flux. In our paper we only deal with the flrst two classes. 

For weak magnetic flux the spectral properties of the system deviate from those predicted 
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by either the GOE or the GUEj^i. In such cases the system is said to be in a crossover^. For 
these systems the Hamiltonian can be decomposed into real symmetric and real antisym- 
metric matrices: 



H 



Hs + iXHa 



(1) 



where X is the crossover parameter— which is equal, up to factors of order unity, to $/$05 
where $ is the magnetic flux through the dot, and = h/e is the quantum unit of magnetic 
flux. Note that the gaussian orthogonal and unitary ensembles are limiting cases of X — > 
and X ^ 1 respectively. 

To understand the meaning of the crossover parameter consider the Aharonov-Bohm 
phase shift picked up by a ballistic electron in a single orbit in the dot: 



For one turn the flux enclosed by the trajectory is proportional to $ = BL^, where L is the 
size of the dot. After N turns the total flux is ^totai = VN^, where factor originates 
from the fact that electron has equal probability to make clockwise or counterclockwise 
orbits, and thus does a random walk in the total flux enclosed. The minimal phase shift 
for the electron to notice the presence of the magnetic flux is of the order 2tt, and thus 
the minimal cumulative flux enclosed by the orbit should be $o = V^^- This leads to 
N = ($o/*&)^; while the time to make N turns is r = LN/vf (for a ballistic/chaotic dot). 
From the Heisenberg uncertainty principle the associated energy scale is: 



where Et is the ballistic Thouless energy^. For a diffusive dot it should be substituted by 
the diffusive Thouless energy = hD/L"^. One can see that when $ is equal to $0) Ex is 
equal to Et which means that energy levels are fully crossed over. 

In this paper the reader will encounter many crossover parameters, and thus many 
crossover energy scales. By a line of argument similar to that leading to Eq. ([3]), it 
can be shown that to every crossover parameter Xi there is a corresponding energy scale 



Ad) = 27r — . 

$0 



(2) 




(3) 
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Breaking the time reversal symmetry of system changes the two-particle Green's function. 
While the two-particle Green's function can in general depend separately on Et, Ex, and 
the measurement frequency u, it turns out that in the universal limit u, Ex <^ E^, it 
becomes a universal scaling function of the ratio Ex/uj- The scaling function describes the 
modification of {G^{E + uj)G^{E)) as one moves away from the "critical" point = 0. The 
limits of the scaling function can be understood as follows: If the measurement frequency 
oj is large (small) compared to the crossover energy scale Ex, the {G^{E + lj)G'^{E)) takes 
the form of the GOE (GUE) ensemble correlation function. If ~ Ex, the Green's function 
describes the system in crossover regime. 

The one-particle Green's function {G^'{E)) is not critical a.s u 0, although it gets 
modified by the interdot coupling. The two-particle Green's function {G^{E + uj)G^{E)) 
always has a diffuson mode^"^, that diverges for small uj in our large-A^ approximation, 
which means that our results are valid on scales much larger than mean level spacing. This 
divergence is not physical and will be cut off by vanishing level correlations for ^ 5 
in a more exact calculation^. On the other hand, the energy scale uj should be smaller 
than Thouless energy of the system for RMT to be applicable. These limitations hold for 
the crossover energy Ex as well. In what follows we study the regime corresponding to 
5 <^uj,Ex < Et. 

The other term that appears in the two-particle Green's function is a cooperon mode. In 
general the cooperon term is gapped if at least one of the crossover parameters is different 
from zero. In the case when the total Hamiltonian of the system is time reversal invariant, all 
the crossover parameters are zero and the cooperon, just like the diffuson, becomes gapless. 
Finally, when each part of compound system belongs to the GUE (the case when all crossover 
parameters are much larger than uj) the cooperon term disappears. 

Our study has a two-fold motivation. The first part comes from works on coupled 
structures with noninteracting particles in acoustic and electronic system3^i^>^, and 
crossovers^i^i^Ui^ii^I. We focus on a complete description of the crossover regimes in all 
three regions (the two dots and the bridge) and define scaling functions for the diffuson and 
cooperon parts of the two-particle Green's function. Using parameters analogous to Ex we 
describe crossover regimes in dots 1 and 2 and the effects of the tunable hopping between 
them. Varying these parameters allows us to obtain results for various physical realizations, 
when different parts of the compound system behave as pure GOE, GUE, or belong to the 
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crossover ensemble. In electronic systems it is easy to break time-reversal by turning on an 
external orbital magnetic flux. In acoustic systems one can break time-reversal by rotating 
the system or a part thereof. As mentioned before, the system of two dots coupled by hop- 
ping has been investigated before using supersymmetry methods^^. However, the authors 
considered only the GUE, whereas here we are interested in the full crossover. In fact, the 
crossover is essential to the second aspect of our work, as will become clear immediately. 

The second part of our motivation is the possibility of using the information gained in 
noninteracting systems to predict the behavior of interacting systemsS^i^iiS^iSS. We consider 
interacting systems controlled by the Universal Hamiltonian2Ii2^i22ii^, which is known to be 
the interacting low-energy effective theory^i^^i^ deep within the Thouless band — ei;'! -C 
Et in the renormalization group^i^ sense for weak-coupling when the kinetic energy is 
described by RMT and the Thouless number g = Et/6 ^ 1. For the GOE the Universal 
Hamiltonian Hu has the forn>2Ii2^i2£i22 



where is the total particle number, S is the total spin, and T = '^Pd'^P,^ ■ ■'■^ addition 
to the charging energy, Hu has a Stoner exchange energy J and a reduced superconducting 
coupling A. This last term is absent in the GUE, while the exchange term disappears in the 



In this paper we concentrate on the reduced Bardeen-Cooper-Schrieffer (BCS) coupling 
A which leads to a mean- field superconducting state when A < 0. Previous work by one 
of us^^ sets the context for our investigation. We consider an interacting system which has 
a single-particle symmetry and a quantum phase transition in the limit Et/S —>■ oo. An 
example relevant to us is a superconducting nanoparticle originally in the GOE. It has the 
reduced BCS interaction and time-reversal symmetry, and the (mean-field) quantum phase 
transition is between the normal and superconducting states and occurs at A = 0. Now con- 
sider the situation when the symmetry is softly broken, so that the single-particle dynamics 
is described by a crossover RMT ensemble. It can be shown^- that this step allows us to 
tune into the many-body quantum critical regime^'^'^ of the interacting system. Thus, the 
scaling functions of the noninteracting crossover are transmuted into scaling functions of 
the interacting system in the many-body quantum critical regime. In our example, the or- 
bital magnetic flux breaks the time-reversal symmetry which is crucial to superconductivity. 





GSE. 
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FIG. 1: The system of two vertically coupled quantum dots. 

When the orbital flux increases to a critical value, it destroys the mean-field superconduct- 
ing state. Above the critical field, or more generically above the critical temperature, the 
system is in the quantum critical regime. 

To be more specific, we consider two vertically coupled quantum dots, the first of which 
has an attractive reduced BCS coupling, while the second has no BCS coupling. Fig. [1] 
shows the geometry, the reason for which will become clear soon. We apply an orbital 
magnetic fiux only through (a part of) the second dot, and observe the effect on the coupled 
system. Our main results are for the mean-field critical temperature Tc of the system, and 
its magnetization in the normal state (above Tc) as a function of the fiux in the normal 
nanoparticle. Such a system could be realized physically without too much difficulty, by, 
for example, growing a thin film of normal metal (such as Au) on an insulating substrate, 
then a layer of insulator which could serve as the hopping bridge, and finally a thin film of 
superconductor (such as Al, which has a mean-field superconducting transition temperature 
of around 2.6K). The orbital fiux can be applied selectively to the Au layer as shown in 
Fig. [1] by a close pair of oppositely oriented current carrying wires close to the Au quantum 
dot, but far from the Al quantum dot. 

The reason for this geometry is that we want to disregard interdot charging effects entirely 
and concentrate on the BCS coupling. The Hamiltonian for the coupled interacting system 
contains charging energies for the two dots and an interdot Coulomb interaction^i. 

+ + U,,N,N, (5) 

Defining the total number of particles as = A^^i + A'"2, and the difference in the number as 
n = Ni — N2 the interaction can also be written as 

re ^ + 16 — + ^^^^ (6) 
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We see that there is an energy cost to transfer an electron from one dot to the other. This 
interaction is irrelevant in the RG sense^'^, but vanishes only asymptotically deep within an 
energy scale defined by the hopping. Our geometry is chosen so as to make Ui = U2 = U12 as 
nearly as possible, which can be achieved by making the dots the same thickness and area, 
and by making sure that their vertical separation is much smaller than their lateral linear 
size. In this case, since is constant, we can ignore charging effects entirely. Charging effects 
and charge quantization in finite systems can be taken into account using the formalism 
developed by Kamenev and Gefen^E^ and futher elaborated by Efetov and co-workers^i^. 
Since our primary goal is to investigate quantum critical effects associated with the BCS 
pairing interaction, we will assume the abovementioned geometry and ignore charging effects 
in what follows. 

After including the effect of the BCS interaction, we find the surprising result that in 
certain regimes of interparticle hopping strength, the mean-field transition temperature of 
the system can increase as the flux through the second quantum dot increases. Indeed, its 
behavior can be monotonic increasing, monotonic decreasing, or nonmonotonic as the flux 
is increased. We can qualitatively understand these effects by the following considerations. 
In the absence of orbital flux, hopping between the dots reduces since it "dilutes" the 
effect of the attractive BCS coupling present only in the first dot. The application of an 
orbital flux through the second dot has two effects: (i) To raise the energy of Cooper pairs 
there, thus tending to localize the pairs in the first dot and raise the Tc. (ii) To cause time- 
reversal breaking in the first dot, and reduce T^. The nonmonotonicity of arises from the 
competition between these two effects. 

Another quantity of interest above the mean-field is the fluctuation magnetization^, 
which corresponds to gapped superconducting pairs forming and responding to the external 
orbital flux. In contrast to the case of a single quantum dot subjected to an orbital flux, 
we find that the fluctuation magnetization^ can be either diamagnetic (the usual case) or 
paramagnetic. A paramagnetic magnetization results from a free energy which decreases as 
the flux increases. The origin of this effect is the interplay between the localizing effect of 
high temperature or the orbital flux in the second dot on the one hand, and the reduced 
BCS interaction on the other. 

The regimes we describe should be distinguished from other superconducting single- 
particle RMT ensembles discovered in the past decade^"^, which apply to a normal meso- 
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scopic system in contact with two superconductors with a phase difference of vr between 
their order parameters^ (so that there is no gap in the mesoscopic system despite Andreev 
reffection), or to a mesoscopic (i-wave superconducting system^. In our case, the symmetry 
of the superconducting interaction is s-wave. However, the most important difference is that 
we focus on quantum critical fluctuations, which are inherently many-body, while the RMT 
classes described previously are single-particle ensembles^>^. 

This paper is organized as follows. In Sec. [TTl we review the basic steps of calculating 
the one-particle and two-particle Green's functions for a single dot. Then in Sec. IIIII we 
present the system of Dyson equations for the one-particle Green's function in the case of 
two coupled dots and solve it in the limit of weak coupling. In addition, we set up and solve 
the system of four Bethe-Salpeter equations for the two-particle Green's function. In Sec. 
|IV]we apply our results to the system of superconducting quantum dot weakly coupled to 
other quantum dot made from a normal metal. We end with our conclusions, some caveats, 
and future directions in Sec. |Vl 



II. REVIEW OF RESULTS FOR A SINGLE DOT. 



Our goal in this section is to calculate the statistics of one and two-particle Green's 
functions for an uncoupled dot in a GOE— > GUE crossover (see appendix |Al and^^ for more 
details), starting from the series expansion of Green's function: 

iP\G-iM.) ^ G^E) ^ (^)^^ . |£ , ^ , ^ , . . . . ,T) 

We are interested in averaging this expansion over the appropriate random matrix ensemble. 
The corresponding Dyson equation for averaged Green's function is: 



> = * + 




The bold line denotes the averaged propagator {G^{E)) and regular solid line defines the 
bare propagator with = E + irj, where 1] is infinitely small positive number. Here 

E stands for self-energy and is a sum of all topologically different diagrams. 

One can solve Dyson equation approximating self-energy only by first leading term and 
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find: 



E\ (9) 

where S is the mean level spacing. This approximation works only for E ^ S. As E gets 
comparable with 5, other terms in expansion for E should be taken into account. 
Then, the average of the one-particle Green's function is given by: 



E 
2 



(10) 



Next, we repeat the procedure for the averaged two-particle Green's function, which can be 



represented by the series: 




+ 



< * c 



+ 



> » > f > 

I I 
I I 

< * < * < 



\ / 

+ X 

/ \ 

< i f < ^ < 



+ 



(11) 

where two bold fines on the left hand side denote {G^{E + u;)G^{E)). The leading contri- 
bution comes from ladder and maximally crossed diagrams. The sum of these diagrams can 
be conveniently represented by Bethe-Salpeter equation. For example, the contribution of 
all the ladder diagrams can be expressed in closed form by: 



n° 



I 

4 



+ 



(12) 



where is a self-energy. For maximally crossed diagrams we have similar equation: 



f 

I 
I 

4 



+ 



(13) 



where 11^ and 11'^ are related to the connected part of two-particle Green's function as: 




+ 



n' 



(14) 



In the limit of cu being much smaller than bandwidth {cu <S N5), the two-particle Green's 
function (connected part) is expressed as: 



{G^JE + u;)GfJE)) ^ 



27r 6apS^s 27r SasS-y/s 1 



+ 



(15) 



The second term is a contribution of maximally crossed diag rams. Ex is a crossover energy 
scale, connected to the crossover parameter as Ex = AX'^NS/n. 

Depending on values of Ex one can speak of different types of averaging. If Ex <^ cu, we 
get average over GOE ensemble, if Ex is of order u, averaging is performed over ensemble 
being in crossover, and, if Ex ^ uj, contribution of maximally crossed diagrams can be 
disregarded, thus going to the hmit of the GUE ensemble. 



III. TWO COUPLED DOTS. 



Next we discuss general framework of our calculation and calculate correlation functions 
for our system of interest, which is two weakly coupled quantum dots (see appendix [Bl for 
more technical details). The Hamiltonian for this system can be represented as: 



ff- = ' : H- 1 : = : ■ d^) 





where Hi and H2 are the Hamiltonians of uncoupled dots 1 and 2. The coupling is realized by 
a matrix V. The elements of Hi, H2, and V are statistically independent random variables. 
We assume that both dots and the hopping bridge are in crossover regimes, characterized 
by parameters Xi, X2, and T respectively. 

In the crossover matrices Hi and V are given by: 

where H^'^ is a symmetric (antisymmetric) part of Hi, and V^'^ is real (imaginary) matrix. 
In what follows we assume that the bandwidths in dot 1 and dot 2 are the same. That 
is, A'^i^i = A^2^2- This should not make any difference in the universal limit 00. 
In addition we introduce the parameter ^ - the ratio of mean level spacing in two dots: 
^ = 61/62- For each realization of matrix elements of the Hamiltonian Htot, the Green's 
function of this system can be computed as follows: 




G = {I®E-H)-'=\ ' ^ ^.z ^^g^ 







) H 


' Gil G12 \ 




[G21 G22 1 



Each element of G has the meaning of a specific Green's function. For example, Gn and 
G22 are the Green's functions that describe particle propagation in dots 1 and 2 respectively. 
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On the other hand, Gu and G21 are the Green's functions representing travel from one dot 
to another. 

Calculating {I ^ E — H)^^ one finds the components of G. For example, 

Gn = [{E - Hi) - V{E - i/a)" V^] = Gi + GiVG2V^Gi + GiVG2V^GiVG2V^Gi + ... 

(19) 

where Gi and G2 are bare propagators in dot 1 and dot 2 defined by Gi = {E — Hi)^^ and 
G2 = {E — H2)~^. 

To find the ensemble average of Gn one needs to average the whole expansion f|T9l) term 
by term. For coupled dots Gij interrelated and in large N approximation can be found from 
the following system of equations: 




The bold straight and wavy lines with arrows represent averaged Green's functions 
{Ga/3,i{E)) and {Gij^2{E)) respectively, while regular solid lines are bare propagators in 
dots 1 and 2. The dotted line describes pairing between hopping matrix elements V, and 
the dashed (wavy) line denotes pairing between matrix elements of Hi {H2). 

The system ( 120|) accounts for all possible diagrams without line crossing. Diagrams 
containing crossed lines of any type are higher order in and can be neglected when 
00. If the hopping between dots is zero, this system decouples into two separate 
Dyson equations for each dot. In the case of weak coupling {U <^ 1), where f/ is a parameter 
controlling the strength of coupling between dots, this system can be readily solved. As zero 
approximation, we use results for a single dot. 

In this approximation one-particle Green's function for dot 1 and dot 2 are calculated as 
follows: 
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E-2T, 







(^) [e + ^VT^ 



1 + ^(1+-^ 



{GUE)) 



(21) 



1 - 



U So 
V? E-2Eo 



(^) [e + ^VT^ 



where e is a dimensionless energy e = -nE /2N5. We used subindex in Sq and {Gq{E)) to 
denote solutions for one uncoupled dot. 

In the large N approximation the contribution to the two-particle Green's function comes 
from ladder diagrams and maximally crossed diagrams. It is convenient to sum them sepa- 
rately. The ladder diagram contribution can be found from the following system of equations: 



n: 



I 

■4- 



+ 



I 

-i- 



n: 



+ 



n: 





where 11^ with proper external lines denote various two-particle Green's functions. As in 
the case of the one-particle Green's function equations, if the inter-dot coupling is zero, the 
system reduces to two Bethe-Salpeter equations for uncoupled dots. 

The system of four equations (l22i) can be broken into two systems of two equations to 
get: 

Dl 
D2 



277 badO- 



-lUO 



9di 



9D2, 
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where go are the scahng functions of diffusion terms in dot 1 and dot 2 defined by: 



1 + 



9di = 



gD2 



(23) 



Here Eu — 2UN5/t^ is the interdot couphng energy scale. These dimensionless functions 
show how diffusion part is modified due to the couphng to another dot. 

Next, for the maximally crossed diagrams the system of equations we have: 



-4- 



+ 





n: 








n: 








n: 





+ 



+ 
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-4- 



n:; 



+ 



+ 



The subsequent solution of this system produces: 
{G^,,AE-ru)Gjp^,{E))c,^ 



27r 
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NiS2 


■ 9 
—lUJ 



9ci 



9C2, 



1 
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1-^— 
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n: 



where gc are the scaling functions for cooperon term defined according to: 



(24) 



(25) 



9c\ 



9C2 



1 I • Ex, +Exo 
1 + I 1 ^ — 



^72 



Ex^Eu ■v/5i?X2-E(7 

1 + 



+ i 



1 I • Ex^ +Ex.2 
1 + I 



Exi Ex2 



Ex^Eij 

~7W 



(26) 
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Here -E'xi,2 = 4X^2^'^/^5 Er = 4r'^Eu/{y^ + are the crossover energy scales, 
describing transition from GOE to GUE ensemble in dot 1 and dot 2, as well as in hopping 
bridge V. 

As we determined how the scaling function gc modifies cooperon part of two-particle 
Green's function and depends on the crossover energy scales defined above, we are ready to 
proceed with write up the connected part of the total two-particle Green's function, which 
is a sum of diffuson and cooperon parts: 



f 

Nf5i -iuJ 1 I - ^xi+gxa _ Ex^Ex2 _ Ex^Eu _ ^/^Ex^Ey / j_\ Eu . ^Er\ 

(27) 

JT^M, l+^y/^^+Z^ ^^^^ 

In general, the coupling between dots changes the bandwidth of each dot. Corrections to 
the bandwidth are of the order of \J and can be neglected for weak coupling. Calculating 
approximations to the second order in V one can be ensure that one-particle and two-particle 
Green's functions can be treated perturbatively. 

Diagrams on Figj2] show the typical behavior of absolute value and phase of scaling 
functions gr, and gc in dot 1. All energy parameters are measured in units of Ey. 

Next we analyze the temporal behavior of the computed statistical characteristics. The 
Fourier transform of the two-particle Green's function shows the time evolution of the density 
matrix of the system. One can observe that the diffuson part of {G^G^) diverges for small 
UJ. To get the correct behavior we replace 1/uj with uj/{uj'^ + rf), and take rj to zero in the 
final result. As for the cooperon term, it stays regular in the small uj limit if at least one of 
the crossover parameters differs from zero. 

First of all, we look at the Fourier transform of {G^G^) in the first dot. We have 
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FIG. 2: Absolute value and phase of diffuson 
1. Frequency lo is measured in units of Eu. 
Ex^/Eu = EX2/EU = 1, Ey/Eu = 0.8, i = l. 



(a,b) and cooperon (c,d) scaling functions in dot 
For these graphs the crossover parameters are: 



+ [2 



iVi 



1 + Yi (g-ia_ _ ^-ta. 



a+ — a. 



(29) 



where a± depend on the crossover parameters (see Eq. fIClip in appendix ICl) 

Then, for the corresponding quantity in the second dot the Fourier transform produces: 
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1 + le-(^+7?)^^* 



a+ — a_ 



• (30) 



IV. TWO COUPLED METALLIC QUANTUM DOTS 



In this section we apply the results obtained in the previous sections to an interacting 
system. We consider two vertically coupled metallic quantum dots, as shown in Fig. [T|, the 
first of which is superconducting and the second noninteracting. For simplicity the quantum 
dots are assumed to have the same level spacing {C, = 1). The calculations presented 
in this section can be extended to the case ^ 7^ 1 in a straightforward way. The first 
(superconducting) quantum dot and the hopping bridge belong to the GOE ensemble. A 
nonzero orbital magnetic flux penetrating the second (noninteracting) quantum dot drives 
it into the GOE to GUE crossover described by the crossover energy scale Ex2- The other 
crossover energy scale Eu describes the hopping between the quantum dots. Because of this 
hopping one can observe a nonzero magnetization in the first particle caused by a magnetic 
flux through the second particle. Roughly speaking, when the electrons in the first dot travel 
to the second and return they bring back information about the orbital flux. 

We wish to compute the magnetization as a function of orbital flux, as well as the mean- 
field critical temperature. It should be noted that since the quantum dot is a finite system, 
there cannot be any true spontaneous symmetry breaking. However, when the mean-field 
superconducting gap Abcs ^ ^) the mean-field description is a very good one^i^i^. Recent 
numerical calculations have investigated the regime Abcs — ^ where quantum fluctuations 
are strong^. We will focus on the quantum critical regime of the system above the mean-field 
critical temperature/field, so we do not have to worry about symmetry-breaking. 

We start with BCS crossover Hamiltonian for the double-dot system including the inter- 
actions in the first dot and the hopping between the dots^^: 
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Hbcsxz — Yl Hfj,^uoc\^s(^uos - XT'^T + 

E eM4,.c^.s - 5\T^T, (31) 

where if*^^-* contains the effect of the orbital flux through the second quantum dot. Here 
T, T"^ are the operators which appear in the Universal Hamiltonian, and are most simply 
expressed in terms of electron creation/annihilation operators in the original GOE basis of 
the first dot (which we call /xq, i^o) as 

Now we need to express the operators c^^^g in terms of the eigenoperators of the combined 
single-particle Hamiltonian of the system of two coupled dots. The result is 

T = Y^ M^^c^,^c^,t5 Mf,^ = V^^(/io)V^^(/io), (33) 

where denotes the eigenvalues of the total system, c^^s operator annihilates electron in the 
orbital state /i with spin s, ip^i^o) is the eigenvector of the compound system, 6 is the mean 
level spacing of a single isolated dot, A > is the attractive dimensionless BCS coupling 
valid in region of width 2ujij around the Fermi energy. Note that while the indices fi, v 
enumerate the states of the total system, the index /iq goes only over the states of the first 
dot, since the superconducting interaction is present only in the first dot. 

To study the magnetization of the first quantum dot in the crossover we follow previous 
work by one of us^^: We start with the partition function Z = Tr{exp—(3H) where (3 = 1/T 
is the inverse temperature. We convert the partition function into an imaginary time path 
integral and use the Hubbard-Stratanovich identity to decompose the interaction, leading 
to the imaginary time Lagrangian 

I 1 2 

= TT - 5^ - + af + aT (34) 

OA 

where a, a are the bosonic Hubbard-Stratanovich fields representing the BCS order parame- 
ter and c, c are Grassman fields representing fermions. The fermions are integrated out, and 
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as long as the system does not have a mean-field BCS gap, the resulting action for o", cr can 
be expanded to second order to obtain 

Seff ^ - fnif3,Ex,uJn)) (35) 

n 

Uf3,Ex,un)= SE\M,A' '-''fX:Zf'^ (36) 

where Un = 2Trn/P, and the sums are restricted to |e^|,|ey| < Hud- We see that the 
correlations between different states fi, v play an important role. Deep in the crossover (for 
Ex ^ 5) we can replace |M^,^p by its ensemble average^^. We will also henceforth replace 
the summations over energy eigenstates by energy integrations with the appropriate cutoffs. 
In previous work— the statistic32^i2ii2^ of |M^i,p was used to obtain analytical results for 
this expression. 

The (interacting part of the) free energy of the system in the quantum critical regime is 
given by^: 

(5F = Eln(l - \f{tUn,P,Ex,)), (37) 

n 

where / is the scaling function given by expression: 

n^u;^,P,Ex,)=5J2\M,.\'^-^^^^f^^^, (38) 

^ + Cj, - tUn 

nu{P) = (1 + exp(/9e,y))~^ is the Fermi-Dirac distribution. We have shifted the energy so 
that the chemical potential is 0. 

Converting this double sum into integral and substituting {M^^]"^ by its ensemble average 
(see Appendices ID] and [E|) . we get: 



rEuT^ (61 - 62)^ + Ex.Eu + El^ tanh(^) + tanh( ^ , 

^" n J-^, ((ei - e^y - Ex.Euf + {Ex, + 2Eune, - e^f ei + €2 - lu^ ' 

(39) 

where ujd is the Debye frequency, and (3 = l/ksT is the inverse temperature. 

One can decompose the ratio in the first part of integrand into two Lorentzians to get^^: 
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fn 



2Ei 



El - El 



In 



+ 



C'/f3^ + {E^ + \uJn\Y 

2 ^2 ^ ^. 



Eu - E^^ 



2Eo 



El 



El 



In 



C'/l3^ + {E2+\uJn\f 



Here C ~ 3.08 and E12 depend on crossover energy scales as follows: 



Elo 



[Ex, + 2Euy - 2EuEx, T J {Ex, + '^Euf{El^ + 4E, 



The magnetization can then be obtained from the free energy: 



(40) 



(41) 



dF 



AL2 dEx, 



9fn 



dB (3 d0 ^i-A/„ 



dE 



X2 



(42) 



where Mnonmt is the contribution from noninteracting electrons^. We will be interested in 
the second term, which is the fluctuation magnetization^. 

For illustrative purposes, we use the parameters for Al in all our numerical calculations, 
with ujd = SAmeV and A = 0.193. This leads to a mean-field transition temperature 
Tco = 0.218me\^ = 2.QK for an isolated Al quantum dot in the absence of magnetic flux. 
In all our calculations we evaluate Matsubara sums with a cutoff exp — \uJn\/ujD- We have 
verified that changing the cutoff does not qualitatively affect our results, but only produces 
small numerical changes. 

It will be informative to compare the two-dot system with a single dot subject to an 
orbital magnetic flux— (see Fig. [3l). We draw the reader's attention to two important 
features. Firstly, the critical temperature Tc decreases monotonically with Ex, resulting from 
the fact that time-reversal breaking disfavors superconductivity. Secondly, the fluctuation 
magnetization is always negative, or diamagnetic, resulting from the fact that the free energy 
monotonically increases as the orbital flux increases. 

Now let us turn to our system of two quantum dots coupled by hopping. Before we carry 
out a detailed analysis, it is illuminating to inspect the behavior of Ei^2 and the coefficients 
of the two logarithms in Eq. (HOl) (which we call A1.2) as a function of Ex,- This is shown 
in Fig. m El tends to Ex,/'^ for Ex, ^ Ejj, and to Ejj in the opposite limit Ex, ^ Ejj. 
E2 tends to Eu for Ex, <^ Eu, while in the opposite limit Ex, ^ Eu E2 — > Ex,- Both 
coefficients 2 start at ^ for small Ex,- For Ex, ^ Eu Ai 1, while A2 —* 0. 
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FIG. 3: Magnetization (per unit volume) in a single dot system as a function of temperature 
for different values of crossover parameters Ex- Panel (d) shows the dependence of the critical 
temperature on Ex- 

The asymptotic regimes T, Ex2 <^ Eu and T, Ex2 ^ Ejj can be understood simply. In 
the first regime, Eu is the largest energy scale, and far below it the spatial information that 
there are two distinct quantum dots is lost. The system behaves like a single large dot with 
a smaller "diluted" superconducting couphng. On the other hand, when T, Ex^ ^ Eu, A2 
is vanishingly small, and the system resembles the isolated first dot with a superconducting 
couphng A but with a crossover energy Eu. Note that the approach of the energies to the 
asymptotes is slow, so for a particular value of Eu it may happen that one cannot realistically 
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FIG. 4: The behavior of Log coefficients in Eq. (I40p and Ei,E2 as functions of the ratio Ex2/Eu. 

approach the asymptotic regime without running into either S at the lower end or uj£, at the 
higher end. Finally, one can envisage situations in which ^ Eu but T > Ejj, for which 
there are no simple pictures. 

The temperature dependence of magnetization per unit volume for different values of 
crossover parameters and Ejj (excluding the part due to noninteracting electrons) is 
shown in Fig. [5l 

In the range where magnetization changes significantly, the fluctuation magnetization 
shows both diamagnetic and paramagnetic behavior. This is in contrast to the case of a 
single superconducting quantum dot subjected to an orbital flux where the fluctuation mag- 
netization is always diamagnetic (Fig. [3l). Close to T = an increase in temperature makes 
the fluctuation magnetization more diamagnetic. A further temperature increase changes 
the fluctuation magnetization from diamagnetic to paramagnetic. For large values of temper- 
ature the fluctuation magnetization is paramagnetic and decreasing as T increases. Another 
set of diagrams. Fig. [6l demonstrates the dependence of the fluctuation magnetization in 
the first dot on crossover parameter in the second dot. Generically, we find that at low 
T the fluctuation magnetization is diamagnetic while at high T it is paramagnetic. 

The variation of crossover energy scales Ex2 and Eu does not change the qualitative 
behavior of the fluctuation magnetization as a function of T or Ex2- A paramagnetic 
magnetization is counterintuitive in superconducting system, because one believes that "an 
orbital flux is the enemy of superconductivity", and therefore that the free energy must 
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FIG. 5: Magnetization (per unit volume) as a function of temperature for different values of 
crossover parameters s^nd Ejj. The fluctuation magnetization is diamagnetic for low T and 
paramagnetic for high T. 

always increase as the orbital flux increases. This assumption is false for our system. The 
explanation is fairly simple, as we will see immediately after the results for have been 
presented. 

The mean-field critical temperature Tc of transition between normal and superconducting 
state strongly depends on and Eu. As one can see from Fig. [3, for very strong 
hopping {Eu ^ Tco) between quantum dots Tc is monotonically decreasing as Ex2 increases. 
On the other hand, for intermediate hopping has a maximum as a function of orbital 
flux, which means that for small values of orbital magnetic flux increases as the orbital 
flux increases. Finally, when Eu is very weak, Tc monotonically increases as a function of 
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FIG. 6: Fluctuation magnetization in the first dot vs crossover parameter in the second dot 
for different values of temperature. The fluctuation magnetization is diamagnetic for low T and 
paramagnetic for high T. 

orbital flux through the second quantum dot. This is in contrast to the behavior of a single 
superconducting quantum dot for which Tc decreases monotonically as a function of orbital 
flux. 

These counterintuitive phenomena can be understood in terms of the following cartoon 
picture. One can think of the two dots as two sites, each capable of containing a large 
number of bosons (the fluctuating pairs). The BCS pairing interaction occurs only on the 
flrst site. When there is no magnetic flux, hopping delocalizes the bosons between the two 
sites, leading to a "dilution" of the BCS attraction and a low critical temperature. The effect 
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FIG. 7: Critical temperature as a function of for several intermediate to strong values (com- 
pared to Tco) of the hopping parameter Eu. For larger values of Ex^ (not shown on graphs) critical 
temperature is equal to zero. 

of the magnetic flux on the second dot is twofold: (i) Firstly, it gaps the cooperon of the 
second dot, which we think of as raising the energy for the bosons to be in the second dot. 
(ii) Secondly, by virtue of the interdot hopping, a small time-reversal symmetry breaking 
is produced in the flrst dot, thereby raising the energy of the bosons there as well. As the 

flux through the second dot rises, the bosons prefer to be in the flrst dot since they have 
lower energy there. The more localized the cooper pairs are in the flrst dot due to effect (i), 
the more "undiluted" will be the effect of the BCS attraction A, and the more favored will 
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FIG. 8: Behavior of critical temperature Tc as a function of for small to intermediate values 
(compared to Tco) of Eu- 

be the superconducting state. However, effect (ii) produces a time-reversal breaking in the 
first dot, thus disfavoring the superconducting state. These two competing effects lead to 
the varying behaviors of Tc and the fluctuation magnetization versus the orbital flux in the 
second quantum dot. When the hopping between the quantum dots is weak {Eu < Tco), the 
first effect dominates, and Tc increases with Ex^- When the hopping is stronger {Eu — Tco) 
the first effect dominates at small orbital flux, and the second at large orbital flux. Finally, 
at very large hopping {Eu ^ Tco), effect (ii) is always dominant. 

When considering the magnetization one must take into account the temperature as 
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well, so the picture is more complex. The general feature is that effect (i) which tends 
to localize the pairs in the first dot also tends to decrease the interacting free energy of 
the system, which leads to a paramagnetic fluctuation magnetization. Effect (ii), which 
breaks time- reversal in the first dot, increases the free energy of the system and thus leads 
to a diamagnetic fluctuation magnetization. Based on our results we infer that at high 
temperature the coherence of pair hopping is destroyed leading to more localization in the 
first quantum dot. The consequences of high T are thus similar to that of the effect (i): A 
lowering of the interacting free energy and a paramagnetic fluctuation magnetization. 

We can make this picture a bit more quantitative for the behavior of with respect to 
Ex- Consider once more the scaling function of Eq. (1401) . which we reproduce here for the 
reader's convenience 



Eu E-^^ + EjjEx2 — 
fn{Ex2,Eu,T) = ^ ^2 _ ^2 1^ 



, Eu El- E\ - EuEx^ 

H n In 

2^2 El - El 



C'/(3^ + {E^ + \u^\)\ 



(43) 



c'ip^ + {E^ + \ujn\y 

It is straightforward to show that reaches its maximum value for uj^ = 0. The condition 
for Tc is then 

~Xfo{Ex2,Eu,T,) = l (44) 

Let us first set Ex2 = 0. Let us also call the mean-field critical temperature of the isolated 
first dot in the absence of a magnetic flux Tco (recall that for the parameters pertinent to Al, 
Tco = 0.218mel^ = 2.6K). Now there are two possible limits, either Eu ^ Tco or Eu ^ Tcq. 
In the first case we obtain 

E'^ 



Tc{Eu)^Ml-j^ + --- ] (45) 



In the second case, Eu ^ Tco, we obtain 



UEu) ^ T,o^e-V^ (46) 

Note that this can be much smaller than Tco and is an illustration of the "dilution" of the 
BCS attraction due to the second dot mentioned earlier. Of course, there will be a smooth 
crossover between the expressions of Eq. (j45l) and Eq. (j46l) . so that Tc is always smaller 
than TcQ. 
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FIG. 9: The behavior of E't^^ vs Eu for numerical simulation and analytical approximation. 

Now under the assumption Ex2, Tc -C Eu we can solve analytically for Tc to obtain 

7?(E.„E,) . -§ + J^e-<^ii-i-'-'m) (47) 

One can further find the maximum of this expression. It turns out that Eu has to be larger 
than a critical value E^ for there to be a maximum. 

Eu = cj^e^^"!^ (48) 

For our values of the parameters uj£, = 34meV, A = 0.193, we find E^ = 0.2A5meV. The 
position of the maximum can now be estimated asymptotically for Eu > E^ as 



Ei^^l6e-E*('^yin|^ (49) 



Figjni compares the dependence of E"^^ vs Eu in case of numerical simulation and the one 
described by Eq. ( H9l) . For large values of Eu compared to E^ the numerically computed 
curve matches the analytical approximation. 
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V. CONCLUSION AND DISCUSSION 



In writing this paper we began with two objectives. We intended to compute nonin- 
teracting scahng functions in the GOE^GUE crossover in a system of two dots coupled 
by hopping, and to use this information to investigate the properties of an interacting 
systent^^i^iiS^iS^ in the many-body quantum critical regime^i^i^. 

We have considered a system of two coupled quantum dots, each of which could have 
its own time-reversal breaking parameter, coupled by a bridge which could also have time- 
reversal breaking. For each crossover parameter, there is a corresponding crossover energy 
scale, which represents the inverse of the time needed for the electron to "notice" the presence 
of that coupling in the Hamiltonian. We have computed the two-particle Green's functions in 
the coupled system in a large- iV approximation^, valid when all energies of interest are much 
greater than the mean level spacing. This allows us to compute the correlations of products 
of four wavefunctions belonging to two different energy levels (which have been previously 
calculated for a single dot for the pure ensembles by Mirlin using supersymmetry methods^, 
and for the Orthogonal to Unitary crossover by Adam et aM). The two-particle Green's 
function splits naturally into a diffuson part and a cooperon part. Each of these parts can 
be represented as times a scaling function, where uj represents the frequency at which 
the measurement is being performed. For example, when we use the two-particle Green's 
function to find the ensemble average of four wavefunctions belonging to two energies, u is 
the energy difference between the two states. The "scaling" nature of the scaling function 
is represented by the fact that it depends only on the ratio of u to certain crossover energy 
scales. For the diffuson part the crossover energy Eu is controlled solely by the strength of 
the hopping between the two dots, while the scaling function for the cooperon part depends 
sensitively on the time-reversal breaking in all three parts of the system. 

In the second part of the paper, we consider the case when one of the dots has an 
attractive BCS interaction, implying that it would be superconducting in the mean-field 
limit at zero temperature if it were isolated, and the other dot has no electron interactions 
but is penetrated by an orbital magnetic flux. The BCS interaction is one part of the 
Universal Hamiltonian2Ii2^i2£i^^ known to be the correct low-energy effective theory^i^i^ 
in the renormalization group^i^ sense for weak-coupling and deep within the Thouless 
band — <^ Et- In order to eliminate complications arising from the charging energy. 
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we consider a particular geometry with the dots being vertically coupled and very close 
together in the vertical direction, as shown in Fig. [H Our focus is on the quantum critical 
regime^i^i^, achieved by increasing either the temperature or the orbital flux through the 
second dot. The first dot is coupled by spin-conserving hopping to a second dot on which 
the electrons are noninteracting. This coupling always reduces the critical temperature, due 
to the "diluting" effect of the second dot, that is, due to the fact that the electrons can now 
roam over both dots, while only one of them has a BCS attraction. Thus, the mean-field 
critical temperature of the coupled system is always less than that of the isolated single 
superconducting dot Tcq. This part of the phenomenology is intuitively obvious. 

However, when the hopping crossover energy Eu is either weak or of intermediate strength 
compared to Tco, turning on an orbital flux in the second dot can lead to a counterintuitive 
increase in the mean-field critical temperature of the entire system. For very weak hopping, 
the mean-field monotonically increases with orbital flux through the second dot, reaching 
its maximum when the second dot is fully time-reversal broken. For intermediate hopping 
strength, the mean-field Tc initially increases with increasing orbital flux to a maximum. 
Eventually, as the orbital flux, and therefore the crossover energy corresponding to time- 
reversal breaking in the second dot increases, the critical temperature once again decreases. 
For strong hopping Eu ^ Tco, Tc monotonically decreases as a function of the orbital flux 
in the second quantum dot. 

We have obtained the detailed dependence of the fluctuation magnetization in the quan- 
tum critical regime as a function of the dimensionless parameters T/Ex2 and EX2/ Eu- Once 
again, the coupled dot system behaves qualitatively differently from the single dot in having 
a paramagnetic fluctuation magnetization in broad regimes of T, Ex2i and Eu. 

We understand these phenomena qualitatively as the result of two competing effects of 
the flux through the second dot. The first effect is to raise the energy for Cooper pairs 
in the second dot, thereby tending to localize the pairs in the first dot, and thus reducing 
the "diluting" effect of the second dot. This first effect tends to lower the interacting free 
energy (as a function of orbital flux) and raise the critical temperature. The second effect is 
that as the electrons hop into the second dot and return they carry information about time- 
reversal breaking into the first dot, which tends to increase the free energy (as a function of 
orbital flux) decrease the critical temperature. The first effect dominates for weak hopping 
and/or high T, while the second dominates for strong hopping and/or low T. Intermediate 
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regimes are more complex, and display nonmonotonic behavior of Tc and the fluctuation 
magnetization. 

It should be emphasized that the quantum critical regime we focus on is qualitatively 
different from other single-particle random matrix ensembles applicable to a normal meso- 
scopic system which is gapless despite being in contact with one or more superconducting 
regions^i^, either because the two superconductors have a phase difference of vr in their 
order parameters^, or because they are (i-wave gapless superconductors^. The main dif- 
ference is that we investigate and describe an interacting regime, not a single-particle one. 
Without the interactions there would be no fluctuation magnetization. 

Let us consider some of the limitations of our work. The biggest limitation of the nonin- 
teracting part of the work is that we have used the large- approximation, which means that 
we cannot trust our results when the energy scales and/or the frequency of the measurement 
becomes comparable to the mean level spacing. When ^ 5 the wavefunctions and levels 
acquire correlations in the crossover which we have neglected. Another limitation is that we 
have used a particular model for the interdot hopping which is analytically tractable, and is 
modelled by a Gaussian distribution of hopping amplitudes. This might be a realistic model 
in vertically coupled quantum dots, or where the bridge has a large number of channels, but 
will probably fail if the bridge has only a few channels. These limitations could conceivably 
be overcome by using supersymmetric methods^ii^. 

Coming now to the part of our work which deals with interactions, we have restricted 
ourselves to the quantum critical regime of the system, that is, when there is no mean-field 
BCS gap. Of course, a finite system cannot undergo spontaneous symmetry-breaking. How- 
ever, in mean-field, one still finds a static BCS gap. The paradox is resolved by considering 
phase fluctuations of the order parameter which restore the broken symmetry To system- 
atically investigate this issue one needs to analyze the case when the bosonic auxiliary field 
a in the coupled-dot system acquires a mean-field expectation value and quantize its phase 
fluctuations. 

We have also chosen a geometry in which interdot charging effects can be ignored. How- 
ever, most experimental systems with superconducting nanoparticles deal with almost spher- 
ical particles. For two such nanoparticles coupled by hopping, one cannot ignore charging 
effects2ii2£iiSili. We expect these to have a nontrivial effect on the mean-field and fluctu- 
ation magnetization of the combined system. We defer this analysis to future work. 
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There are several other future directions in which this work could be extended. New 
symmetry classes^*^ have been discovered recently for two-dimensional disordered/ballistic- 
chaotic systems subject to spin-orbit coupling^i^. In one of these classes, the spin-orbit 
coupling is unitarily equivalent to an orbital flux acting oppositely^"^ on the two eigen- 
states of a single-particle quantum number algebraically identical to cr^. Due to the unitary 
transformation, this quantum number has no simple interpretation in the original (Orthog- 
onal) basis. However, it is clear that the results of this paper could be applied, mutatis 
mutandis, to two coupled two-dimensional quantum dots subject to spin-orbit couplings. In 
particular, consider the situation where one quantum dot has no spin-orbit coupling, but 
does have a Stoner exchange interaction, while the other dot is noninteracting, but is made 
of a different material and has a strong spin-orbit coupling. Work by one of us has shown^^ 
that by tuning the spin-orbit coupling one can access the quantum critical regime, which is 
dominated by many-body quantum fluctuations. The above configuration offers a way to 
continuously tune the spin-orbit coupling in the first dot by changing the strength of the 
hopping between the dots. 

In general, one can imagine a wide range of circumstances where changing a crossover 
parameter in one (noninteracting) dot allows one to softly and tunably break a symmetry 
in the another (interacting) dot, thereby allowing one access to a quantum critical regime. 
We hope the present work will be useful in exploring such phenomena. 
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APPENDIX A: ONE UNCOUPLED DOT 

In this Appendix we calculate one-particle and two-particle Green's functions for a single 
dot undergoing the crossover. The strength of magnetic field inside the dot is controlled by 
crossover parameter X. The Hamiltonian of the system in crossover is: 
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(Al) 

where Hs^a are symmetric and antisymmetric real random matrices with the same variance 
for matrix elements. Normalization (1 + X^)~^/^ keeps the mean level spacing 5 fixed as 
magnetic field changes inside the dot. 

We define the retarded one-particle Green's function as follows: 



^E+-Hj^p E+ \ ' E+ ' (E+y E+ ' {E+f {E+f 

(A2) 

Here if is a Hamiltonian, and is the energy with infinitely small positive imaginary part 
E+ = E + ir]. 

This series has nice graphical representation: 

I I I 

I I I 

G^{E) = • + — i — * — + — > — '^-^ + .., (A3) 

where straight solid line represents 1/E^ and dashed line stands for Hamiltonian. 

= > , H = 1 (A4) 

Just as in disordered conductor or quantum field theory the target is not the Green's 
function itself, but rather its mean and mean square. We take on random matrix ensemble 
average of Gais- Such averaging assumes knowledge of (if"), where angular brackets stand 
for gaussian ensemble averaging, and n = 1, oo. For = 1 we have (H) = 0, while for n = 2 
the second moment reads: 



All higher moments of H can be computed using Wick's theorent^. Thus, the ensemble 
averaging leaves only the terms containing even moments of H. Introducing the notation 
for {HH) = we obtain, for the averaged series: 



{G^{E)) = > + > > > + . . . (A6) 
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Then, the expansion (1A6|) can be written in a compact form of Dyson equation: 
> = > + > 




(AT) 

The bold hne denotes the full one-particle Green's function averaged over Gaussian en- 
semble, and S is a self-energy, representing the sum of all topologically different diagrams. 
The corresponding algebraic expression for the Dyson equation can be easily extracted from 
Eq. flA7p producing: 

Ga/S = ^ Gau^iyfM^ + (A8) 

UfJ, 

where Gap means (Gf^ (£")). Now, using the fact that Gai3 = GaSap and T^ais = ^a^ap (no 
summation over a implied), one can solve this equation and obtain: 

Ga, = (A9) 
Next we approximate self-energy by the first term in large N approximation: 



7 ^ 



N5V 5, 



J 



(AlO) 



Solving Eq. (1A10I) for the self-energy we determine: 



E i \f2N5V 



Consequently, the ensemble average of one-particle Green's function is given by: 



{G%(E)) = ; {G%(E)) = {G%(E))\ (A12) 

f + i\/(¥)^ 

Next, to study the two-particle Green's function we notice that the main contributions 
come from ladder and maximally crossed diagrams: 




> f > > f > f > > » > j i > 



< < * < < i < * < < i f < ) i < 



(A13) 
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Two bold lines on the left side stand for the average two-particle Green's function {G^{E+ 
ijj)G'^{E)). The sum of ladder diagrams is described by Bethe-Salpeter equation: 
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(A14) 



(A15) 



where is a ladder approximation of diffuson part of two-particle Green's function. Here 
F[E, to] is a product of two inversed averaged one-particle Green's functions and in the limit 
u ^N5 is: 



F[E,u;] = {G\E + uj))-\G^{E))-' 



iu6N (N6\- 



One can solve this equation taking into account H^;^'^ = H^SasSis-y' 



F[E,u;] 



(A16) 



9- (A17) 
F[EM-{^T 

Multiplying 11^ by F'^[E,u] we arrive at the following expression for the diffuson term: 



{GUE + u;)GfJE))n- '^"^'^^^ 



(A18) 



Then, we turn our attention to the equation for maximally crossed diagrams. We have 
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(A19) 



and n*^ is expressed in terms of F[E, lo] again: 
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(A20) 



1 + Xy FiEM-mi'ff 
Assuming X to be small compared to unity (weak crossover) , we evaluate the contribution 
of maximally crossed diagrams to Green's function to get: 



-iuj 1 + 



(A21) 
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where Ex — 4:X'^NS/7r is a crossover energy scale. Final expression for the connected 
part of the two-particle Green's function is: 



APPENDIX B: TWO COUPLED DOTS 

This Appendix contains details of the derivation for statistical properties of the Green's 
functions for the two coupled dots connected to each other via hopping bridge V. Coupling 
between dots is weak and characterized by dimensionless parameter U. For the system of 
uncoupled dots the Hilbert space is a direct sum of spaces for dot 1 and dot 2. Hopping V 
mixes the states from two spaces. The Hamiltonian of the system can be represented as: 




Htot^ I .. I ■ (Bl) 
For Hi^2 and V we have: 

Here S (A) stands for symmetric (antisymmetric), and R (I) means real (imaginary). 
Below we use Greek indices for dot 1, and Latin indices for dot 2. We also found it convenient 
to keep bandwidth of both dots the same; that is, A^i^i = A^2(^2 with ^ = 5i/52- 

The following averaged products of matrix elements of H can be obtained: 



(B3) 



where Xi and X2 are the crossover parameters in dot 1 and 2. Pairings between V matrix 
elements are: 

... . .t X VNMS2U. . 
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with r a crossover parameter in hopping bridge. Normahzation for V pairing is chosen to 
coincide with that of (HH) when = 1. 

To determine one-particle Green's function we use the system hsted in Eq. (1201) . The 
straight and wavy bold lines with arrows represent averaged functions (Gf(£')), (G'f^(-E')) 
in dot 1 and 2, regular lines represent bare propagators, and the rest of the lines describe 
pairings between Htot matrix elements. We have: 



1 

E+ 
1 

~~ E+ 



(B5) 



I 



The corresponding analytical expressions of this system of equations are: 



H H 



1 



^2 ^+ ^ ^+ ^^+: 



with Gi and G2 connected to Green's functions via: {G^^^-j^{E)) = GiSa-y, (G^gl-^)) ~ G26i, 
The self-energies E"™ are to be determined using standard procedure^. 
We observe, that the system of two linear equations ([B]) has a solution: 



Go 



1 



Here we approximated self-energies by the first term in large N expansion again. In this 
approximation evaluation of S""^ yields: 



J a 13 



E+ - - S12 



Vl2 _ ■^12r 
^af3 — ^ "qI3 



^22 



^21 



^f^ = G2j2{v^.v;i,) 



VNj^N26i62U 6, 



af3 



E+ - S22 - S21 



N262 



TT 



E^ 



■^22 



VWN2N16162U 



E+ 



^21 



^11 



n2 ' 
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Thus, to find all S"™ one needs to solve the following system of equations: 











_ J.22 _ 








f/S^Vv^ and S 


12 _ 



x/]Vi7V2A^2'5i52t/ 



(B6) 



VWN2N161S2U 



jj-2 



( ]B6p . For example, the pair of first and third equations can be rewritten as: 



2 



2 



(B7) 



For weak coupling the solution can be found by expanding self-energies S^^ and S^^ in 
series in U. Taking the solution for single dot as zero approximation (below all the solutions 
for the uncoupled dot will be marked with subscript 0) we get 



= J^l^ + UEl^ (B8) 
= s22 + ?7Sf. (B9) 

Note that A^i^i = N262, and Sj^ = Sg^ = Sq. 

Plugging into the right hand side of Eq. flB9P in system flB7l) we arrive at: 



J.21 ^ U^u 
= f/^S^^ 



v^^+-2So; (BIO) 
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Neglecting the higher powers in U for one-particle Green's functions we finally arrive at 
the following expressions for the single particle Green's functions: 



{G%,om 

I _ U So 



1 + 



'VT- 



6ij 



v/|E-2Eo 



(^) [e + ^VT^ 



where e = 7rE/2N5. 

Now we switch our attention to the calculational procedure for the average of the two- 
particle Green's functions {G^^^E + uj)Gff^^^{E)) and {G^^E + u)G^^^^{E)) . In the limit 
of large A'^i and ladder and maximally crossed diagrams contribute the most. For ladder 
diagrams we obtain the system of Bethe-Salpeter equations (see Eq. ( !22l) ). Here we used 
the following notation: 



n: 



{G^,,{E + u)G^,{E)), 



n: 



{GUe + uj)GUE)) 



n:: 



{GUE + uj)GUE)), 



n: 



(Bll) 



{Gl{E + uj)GUE)) 



For the diffuson n£„ the system of algebraic equations reeds: 



" 7r2 ~^ 7r2 Fi[E,uj] tt^ F2[E,uj] 

_ N26I N26I N2Ug , ./Njr26i62U iViHf^ 



2' 7r2 7r2 F2[E,u;] Fi[E,u;] 



^2 7r2 ^ 7r2 7r2 F2[E,a;] 



(B12) 



2' 7r2 ^ 7r2 F2[E,cu] 7r2 Fi[E,cu]' 

where cj] and F2[E, oj] are defined as products of inverse averaged one-particle Green's 
functions in the first and second dots respectively. For small values of U and uo these 
functions can be approximated as follows: 
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,-1 



F,[EM = {GnE + u;)r\Gf{E)) 
F,[E,u;] = {G^iE + u;))-\G^{E)r' 



TT 



N2S2 



TT 



2 r 



1 + v^f/ - iCo 

1 + -^-.. 



(B13) 



where lj = Tru!/2N6. The system of four equations given by the Eq. (]B12p can be decoupled 
into the two systems of two equations each. To determine one solves the system of the 
first and the last equations of Eq. (1B12I) to get: 



Fi{E,u) 



n 



n 



^ ^Njf2Ni6i62U n 



^fWN'25i52U 



(B14) 



F2{E,uo))^^^ 7r2 F^{E,u) 

Then, solving the resulting system (Eq. flB14p ) and attaching external lines one obtains 
expression for the two-particle Green's function in dot 1: 



{Gl^AE + uj)G%AE))o 



27r 5nn5. 



1 J i Eu 



The corresponding correlator for dot 2 is readily obtained as well: 



(B15) 



{Gl2{E + uj)Gt,,2{E))D- ^,2 



27r 5;,5ik 



(B16) 



For the second part of the Green's function (which is the sum of maximally crossed 
diagrams) the system of equations is described by Eq. (!24l) . Transforming this graphical 
system into the algebraic one, we get: 



■'-'■22 



■'-'^12 



n: 



c 

21 



l-Xl\ Ni5f fl-Xf\ Ni5f NiU 



I- XI 



XI 



1 + r2 
i-r^ 
TTT2 



7r^ 
X2SI 
7r2 



■■11 



1 



XI 



TT 

N2Si X2^'i2 



F^E.uj] 



+ 



i + r2 



\-Y'^\ v^iViiVa^i^sf/ iV2n^i 

i-r2 



F2\EM 



c 

12 



TT 



+ 



VWN26162U 



Y2 
^2 



1 + r2 



7r2 F2[E,uj] 
7r2 F,[E,iu] 



+ 



+ 



Fi[E,Lu] 

TTxf 



Fi[E,uj] 



i-xl 

1 + 



N251 iVsH^i 



F2[E,uj] 
(BIT) 
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Once again, the system at hand breaks into systems of two equations each. We proceed 
by combining the first and the last equations to obtain: 



1 - 



n 



l + XlJ F,[E,cu] 
T+p; Fi[E,u;] 



1 - y/Nj7^N2SiS2U 



21 



2 



1-Xf\ Ni6j 



+ 



1 - 



I + XIJ F2[E,uj] 



^''21 



l + XfJ 7r2 

1 - y/Nj^5iS2U 



1 + r2 y TT- 
(B18) 

Now we can construct approximations for the expressions, containing crossover parame- 
ters. For example, for small values of X and F the solution for U^^ is expressed as follows: 



i - 2X2) (-^ 2X1) + (1 - 4F2)f/2 



(V^t/ + 2X2) (-^ _ ^ 2X|) - (1 - 4F2)t/2 
Next, introducing crossover energy scales: 



(B19) 



E, 



AX' 



.NS 



E, 



TT 



U 



2U 



NS 



TT 



we obtain the solution for Yl^i in the following form: 



Er 



AT^Eu 



(B20) 



c _ NiSf 1 



1 - 



E 



_£2 



11 



TT^ — iuj -1 +^X2 _|_ Exi Ex2 _|_ Ex^Eu sfjEx^ Eu 



(B21) 



Then, adding external lines to nf^^ for Green's function we get: 



(Gf,,i(^ + ^)G/^,i(£;))c = 
27r (5o,g(^-y^ 

N^bx —iuj 1 I A ^x-^+Ex^ _ Ex^Ex^ _ Ex^Eu _ -/lEx^Ey . ( 1 A /-. , • \ ' 

(B22) 

Similar manipulations for the corresponding correlator of Green's functions for the second 
room result in: 
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(B23) 

Finally, the connected part of the total two-particle Green's function is obtained as a sum 
of diffuson and cooperon parts, yielding: 



9 r r X _|_ J-Eu. 

N?Si —iuj 1 I ,- gxi+gx2 _ ^^1^X2 _ Ex^Eu _ ^/lEx2Eu ■ ( ij \ \ Eu 1-^ , ■E-c\ 

(B24) 



27r 1 + 



27r Mi,- 



(B25) 

APPENDIX C: FOURIER TRANSFORM OF TWO-PARTICLE GREEN'S 
FUNCTION 

To be able to study temporal behavior of electrons in the rmt system we introduce the 
Fourier transform of two-particle Green's function. We define it via the following integral: 

1 r+oo 

{GUt)G%(t)) = y ^ exp— + uj)G%{E))dujdE. (CI) 

To get the correct behavior of the diffuson part for small u, we replace l/uhy u/ {uj'^-\-rf) , 
where 77 is infinitesimal positive number. Now we introduce for dot 1: 
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2vr dais^-rS 



111 Eu 



27r 



io; ( a; + -^E] 



u 



The Fourier transform of this diffuson term gives: 



(C2) 



foit) 



1 

2^ 



exp{—iu!t) fD{i^)diu ■ 



(C3) 



Next steps are the standard steps of integration in complex plane. For t > one closes 
contour in lowerhalf plane. One root is located in upper half plane and two more are located 
in lower half plane. The integration yields: 



2ti 



As T] approaches zero, fDit) becomes: 



+ 



{Vl+^)VlEle 



iv^+^yE'u-v' 



(C4) 



2% 



2 



'N?s,i+a^'^^ J' ^^^^ 

The full Fourier transformation includes integration over E as well. In current approx- 
imation, when E is close to the center of the band, (Gf G^) is independent of E. It will 
depend on E if we integrate over the whole bandwidth. The exact dependence of (GiGf) 
on E far from the center of the band is not known. To get correct expression we assume 
that integration over E adds to (Gf-Gf) multiplicative factor A^i^i along with normalization 
coefficient A. Also, for index pairing a — (3 and 7 = 5, G^^Gf^ becomes transition proba- 
bility density P(t)a^^- Using equipartition theorem, for t — > oo summation of P{t)a^^ over 
a one can get total probability to stay in dot 1. It is equal to Ni/{Ni + N2). That is. 



iVi 



(C6) 



N1 + N2 1 + e 

Integration over E and summation over a gives the factor of AN^Si. We identify the 
normalization constant as A — I/tt. Note, that we did not use cooperon part fc{t) to 



42 



determine normalization constant A. The reason for that is chosen index pairing. After 
the summation over a cooperon part contribution is of the order 1/Ni compared with the 
diffuson part. After integration over E with proper normahzation foit) becomes: 



2 ^ 



iVi(l + 

Then we perform the Fourier transform of the cooperon part: 



(C7) 



fciu;) 



J_Eu 



N?6i —iuj 1 _ Ex^+Ex^ Ex^Ex^ Ex^Eu VlEx^Ey ( j_\ Eu_ ( Er _ t\' 

(C8) 

The fc{'^) is a regular function when uj approaches limiting values, provided at least one 
of the crossover energy scales Ex^-, Ex2, or E^ differs from zero. 

To make fc{'^) more suitable for the Fourier transform we manipulate Eq. (ICSP into: 
27r 



TP 



i^f - {{Ex, + ^^^Eu) + {Ex, + {lu) 

+ {ex.Ex, + + Ex,Eu^i+ ( + -^)EuE^ 

and observe that the poles of fc{<^) are given by 



(C9) 



{Ex, + V^Eu) + {Ex, + ^) ± 



(CIO) 



with D = {{Ex, + V^Eu) - {Ex, + Eu / Vl)f + 4^^(1 - AT^). The parameter D is always 
positive and uj± are imaginary complex numbers. 

It can be proved that {Ex, + \/iEu) + {Ex, +Eu/ \f^) > \^ for all values of parameters, 
which means that the poles are pure imaginary numbers in lower half complex plane: 



{Ex, + VlEu) + {Ex, + f|) ± 



(Cll) 



The function fc{'^) now reads: 



2n iuJ -Ex,-^ 



N'^di {iuj — a-){iuj — a+) 



(C12) 



43 



We perform the Fourier transform and use the normahzation factor to obtain: 



fc{t) 



27r 



exp{—iLjt)fc{uj)dujdE = SasS- 



7/3 



1 + 



a+ — a. 



(C13) 



Hence, the full expression for the Fourier transform for the two-particle Green's function 
in the dot are given by: 



iVi(l + [2 



1 + ^^-i^+^)Eut 



1 + ^ yi (g-ta. 



a+ — a_ 



(CM) 



22 



2i 



iV2(l + 



1 + l^-iVl+^)Eut 

2 e 



No 



a+ — a_ 



, (C15) 



where a± is defined through Eq. fIClip . 



APPENDIX D: CORRELATION OF FOUR WAVE FUNCTIONS 

In this appendix we obtain correlation of four wave functions {ipn{.Oi)ilj'^{;'y)ipm{^)'^mil^)) 
for the system of two coupled dots. This has been obtained in a single dot for the pure 
ensembles by supersymmetry methods by Mirlin^'', and for the GOE^GUE crossover by 
Adam et at^. We consider ensemble average of the following product: 



{[G%{E + u:)- Gt^{E + u)] [G%{E) - G%{E)\) ^ 

- {G^^{E + uj)G%{E) - {Gi^iE + uj)G%{E)) = -2{5^p5,sRe[D,] + 5^s5,pRe[C,]), (Dl) 



where Di and Ci are the diffuson and cooperon expressions from Eq. fl271) . Here we used 
the fact that ensemble average of G^G^ and G^G^ are smaller than G^G^ and G^G^. 
On the other hand, we have: 
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GU^) - Gt,{E) = -2m J2 Mc^mmE - E„), (D2) 

n 

and 



( [G^^iE + u;)- Gt^{E + uj)] [GI{E) - G%{E)] ) ^ 

- 47r2(^^„(a)V;:(7)V'n^(5)C(/5)5(i^ + a; - En)5{E - E^)). (D3) 

71, m 

We know that in the crossover components of eigenvalues and eigenvectors are correlated 
with each other. This correlation is small already on the distances of a few 5 and can be 
neglected in the limit oj ^ 5, so Eq. fIDSP can be approximated by: 

-A7:^i^n{a)rMi^in{5)rnm{Y.^{E + uj - En)){Y,5{E - E^)). 

n m 

where n and m mark energy levels close to E + u and E respectively. 

The average of the sum is a density of states p{E) = — E„)) = 1/5. Then, we 

get 



For the two coupled dots we have: 

In order to calculate i?e[Ci] from Eq. ([27]) we are going to assume that magnetic field is 
zero in the first dot and in the hopping region {Exi = E-p = 0), and the second dot is in 
GOE to GUE crossover {Ex2 ~ t^)- Then, 

R ir]= 27r V^EuLu' + (Eu + ./^Ex,)EuEx, , . 

NfS,{u;^-V^EuEx,y + {Ex, + {V~^ + j^)Euyu;^ ^ ^ 

The relation between the mean level spacing 6 for the system of coupled dots and the 

mean level spacing in the first uncoupled dot 6i is as follows. The averaged density of states 
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in coupled system is going to be the sum of densities of each dot: (p) = (pi) + (P2), or 
5^^ = + 52^"'^. Thus, we conclude that 5 = 5i/ {I + 

Finally, we set Eq. (IDip and Eq. (lD4p equal and obtain correlation of for the wave 
functions: 

APPENDIX E: SUM RULE FOR DOUBLE DOT SYSTEM 

To verify the expressions we have obtained for the averaged Green's functions we use a 
sum rule. 

The pair annihilation (creation) operator T{T'^) in the basis of two uncoupled dots is a 
sum of two terms belonging to each dot: 

~ ^ ] '^"o.-l'^Qo.T ~^ ^ ^ ^io,l^'i-o,h 

(El) 

Greek indices go over the states in the first dot, and Latin indices go over the states in 
the second dot. The subindex denotes the basis of two uncoupled dots. 

Our first goal is to calculate the commutator [T^,T]. As operators from different dots 
anticommute, one gets: 

[T\ T]=Y, [^ioAo^V ^PoACpoA + E[i,T<a' c.oaC,o,T] = ^le + N,, -N,- N„ (E2) 

ao >/3o «0 JO 

where Nie, N2e are the operators of total number of electrons in dot 1 and dot 2, and A^i, A'2 
are the total number of levels in dot 1 and dot 2. 

The expectation value of [T^,T] in ground state at zero temperature is: 



[T\T] = {n\[T\T]\n) = N^- N. (E3) 



46 



A^e and N are the total number of electrons and levels in both dots. This number is conserved 
when going to another basis. 

Now we choose the basis of the system of coupled dots. In this basis Cag^s — 
'^m('^o)Cm,s, and = '^„^'ipmi'io)Cm,s, where Cms is annihilation operator in new basis. 
Using this transformation, we rewrite pair destruction operator as follows: 

^ — ^ (^ao,lCao,1 + ^ Qo,J.^io,T ~ -Dmim2Cmi,J,C^2,T' (^^) 

ao in mi,m2 

where Dmim2 is defined by the following expression: 

Dmim2 = ^ J^V'mi(Qio)'0m2(Q;o) + ^'0mi(«o)'0m2(«o) Cmi,J.Cm2,T 
mi,m2 \ ao io / 

= Yl bo) V'ma (Po) ■ (E5) 
Po 

The index po runs over all states in the first and second dots for the basis of uncoupled dots. 
In the new basis the T, operators look like this: 



TOim2Cmi,J,Cm2,T 



(E6) 



mim2 ra2,'\ ra\,\,' 

mi, 1712 



Consequently, in the new basis. 



[T^,T] ^mim2-^m3m4[(^ln2,'\^mi,l^^rn3,i^rn4,'\] 

mi ,7712 7773,7774 

~ I -^77717772-^^717774 j '^^S ,T ,T ~ I -^777 1 7772 -^"^3^2 l '^"^3 , J. , J ■ (E7) 

7772,7774 \ 7771 / 7771,7773 \ 7772 / 



One can go further and use completeness condition Xlm '^m(Po)'^777(^o) = <^po77o show 
that in the new basis the value of commutator is N^, — N. Our next goal, however, is to take 
the disorder average of the vacuum expectation value and to prove the invariance of [T^ , T]. 

Taking into account that {^\cl^j^^^Cm2,'[\^) = (^77717772 ©(a« - -£^7771) and {^\cm2,i(iii,i\^) ^ 
(^77717772(1 ~ Q(a* ~ -E'7771))) the ground state expectation value for the commutator is: 
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[Tt,T] = {n\[T\T]\n) = J2 \Dm,m,\'[2Q{fi-E^,) - 1], (E8) 

nil i»^2 

where 9(x) is a step function. 
Averaging over disorder gives: 

mi,m2 mi, 7712 

Converting this into integral, we get: 



/ dE,dE2p{E^)p{E2){\D{E,,E2 
-w J-w 



W rW 

dEidE2p{Ei)p{E2){\D{Ei,E2)\^) (ElO) 

w J-w 



The density of states p{E) is the Winger's semicircle law: 



y /y 



2N 

where 2iy is the bandwidth and is the number of states in the system. 
To proceed we need to find the ensemble average of the following object: 



?) = '^{i^*mi{PoWm2iPo)^nii{no)^m2{no))- (Ell) 



PO,no 



Using results of appendix [D] one can obtain expression for the correlation of four wave 
functions in the form: 



pono 



27r2 

-J2(^noPoiE2)GZ,pSE,))\. (E12) 

Pon-o 

Note, that to get the correct answer for the sum rule one should keep (G^G^) term as 
well. Summation in Eq. (lE12p is performed over the states in both dots. 

When the dots have equal mean level spacing 61 = 62 = Sq, one particle Green's function 



can be found exactly from the system (IB7p without approximation in U : 
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where W = 2No5oV^ + U /vr is the half bandwidth and sin0 = E/W. Here both indices 
Po and p'q belong either to the first or to the second dot. 

The sum in Eq. (1E12|) can be broken into four sums, when the indices Pq, Uq belong either 
to the first dot, or to the second dot, or one of the indices go over the states in the first dot, 
and the other one goes over the states in the second dot. 

For example, for {G^G^) part we have the following expression: 



Pong 



Wj ' ' [(1 + t/)e~'^<^2i _ 1] [(1 + f/)e-i<^2i - C] - 
ATlAVn (1 + U)e-''^^^ - 1 

+2No(^] (1 + f/)- ^ 



(E14) 



Wj ' ' [(1 + f/)e-^-^2i - 1] [(1 + f/)e-'<^2i -C]-U^ 



Here 02i = 02 - (pi, and C = (1 - + ^1) 

The first term in Eq. ( 1E14I) is the contribution of {G^){G^) plus the cooperon part of 
two particle Green's function in the first dot. The second term describes contribution of 
free term and cooperon part in the second dot. The last term is a sum of transition parts 
from dot 1 to dot 2 and vice versa. It appears that these transition terms are equal, which 
explains coefficient 2 in front of the last term in Eq. (]E14p . 

Summation of the (G^G^) gives similar result: 
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Po"o 



WJ ' ' 1(1 + U)e-<*=' + 1]|(1 + U)e-'*" +(]-U^ 



Wj ' ' [(1 + f/)e-*'^2i + 1] [(1 + f/)e-i'A2i + (] - [/2 
2\^, U 



VvrJ ^ [(1 + U)e-'^^^ + 1][(1 + U)e-'^^^ + C] - f/'' 
where ■?/'2i = 02 + 0i- 

In principle, there should be terms corresponding to diffusons in dot 1 and dot 2. However, 
these terms after summation over po,nQ are smaller than the others and in the large 
A^o limit can be neglected. 

Although one can use Eq. flElOp to verify the sum rule, it is more convenient to work 
with derivative of Eq. (lElOp over /i at /x = 0. 

It gives: 

— ([Tt,r])^=o = 2p(0) / dE2p{E2){\D{E, = 0,^2)1'). (E16) 

Cp J-W 

On the other hand, this expression should be equal to: 

— (AT, - AT) = — (2- + 2 p{E)dE -N)= 2p{p). (E17) 

Comparison of Eq. flE16l) and flE17l) at /i = results in the following condition for the 
sum rule: 

rW 

/ dE2p{E2){\D{E, = 0,^2) H = 1. (E18) 
J-w 

The integral in Eq. (lElSP was computed numerically and matched the unity with high 
accuracy. 
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